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Abstract 

Current neutrino oscillation data from solar, atmospheric, and reactor experiments 
are consistent with the neutrino mixing matrix elements taking values sin 2 #12 = 
1/3, sin 2 #23 = 1/2, and sin 2 #13 = 0. We present a class of renormalizable gauge 
models which realize such a geometric mixing pattern naturally. These models, which 
are based on the non-Abelian discrete symmetry A4, place significant restrictions on 
the neutrino mass spectrum, which we analyze. It is shown that baryogenesis via 
leptogenesis occurs quite naturally, with a single phase (determined from neutrino os- 
cillation data) appearing in leptonic asymmetry and in neutrinoless double beta decay. 
Such predicted correlations would provide further tests of this class of models. 



Introduction 



Our understanding of the fundamental properties of neutrinos has improved dramatically 
over the last few years. Atmospheric and solar neutrino experiments have by now firmly 
established occurrences of neutrino flavor oscillations pQ. In order for neutrinos to oscillate, 
they must have non-degenerate masses. In addition, different neutrino flavor states must 
mix with one another. When positive evidence for oscillations from solar and atmospheric 
neutrinos are combined with results from reactor data [21 Ej, one obtains the following 
neutrino mass and mixing pattern (with 2a error bars) 

Am 2 =m\-m\ = 7.92 x l(n 5 (l ± 0.09) eV 2 , (1) 
Am 2 tm = m\ - ml = ±2.4 x 10- 3 (l+°; 2 }) eV 2 , (2) 
sin 2 12 = 0.314(11°;}*) , sin 2 6 23 = 0.44(1^) , sin 2 9 l3 = 0.9±°; 23 x 10~ 2 . (3) 

Here are the (positive) neutrino mass eigenvalues, and Oij are the neutrino mixing angles. 
m\ — m\ > in Eq. (1) is necessary for MSW resonance to occur inside the Sun. The sign 
of Am 2 tm , which is physical, is currently unknown. 

A remarkable feature of the oscillation data is that they are all consistent with a "ge- 
ometric" neutrino mixing pattern defined by the parameters sin 2 #12 = 1/3, sin 2 #23 = 1/2, 
and sin 2 9 13 = 0. In fact, these geometric mixing angles are very close to the central values 
of Eq. (3). We observe that unlike the quark mixing angles, which are related to the quark 



mass ratios in many models (eg: 6c — ymd/m s ), the neutrino mixing angles seem to be 
unrelated to the neutrino mass ratios. 

The purpose of this Letter is to provide a derivation of such a geometric neutrino mixing 
based on renormalizable gauge theories. The neutrino mixing matrix (the MNS matrix) that 
we will derive has the form j3J EJ EJ [7| 



U 



MNS 



(n A \ 

1 1 i_ 

\ L JL 1 I 

X V6 V3 V2 ' 



P. (4) 



Here P is a diagonal phase matrix which is irrelevant for neutrino oscillations, but relevant 
for neutrinoless double beta decay. Eq. (HJ) yields the desired values, \U e 2\ 2 = 1/3, {U^l 2 = 
1/2, \U e3 \ 2 = 0. 
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Our derivation of Eq. (@J) will be based on the non-Abelian discrete symmetry A4, 
the symmetry group of a regular tetrahedron. This symmetry group has found application 
in obtaining maximal atmospheric neutrino mixing \8 ; and in realizing quasi-degenerate 
neutrino mass spectrum [Oj. No successful derivation of Eq. has been achieved to our 
knowledge (based on A4 or other symmetries) in a renormalizable gauge theory context. For 
attempts along this line see Ref. |U [TUl HI]. In Refs. |U Ej, Eq. (JU) was suggested as a 
phenomenological ansatz. In Ref. [10], a higher dimensional set up is used to motivate Eq. 
(jlj). Ref. ^T] analyzes special cases of an A4 derived neutrino mass matrix towards obtaining 
the structure of Eq. (j3J). A large number of models in the literature have derived maximal 
atmospheric mixing based on non-Abelian symmetries ^2]; but in most models the solar 
mixing angle is either maximal (now excluded by data) or is a free parameter. 

We will see that the derivation of Eq. places strong restrictions on the neutrino mass 
pattern. We find that the out of equilibrium decay of the lightest right-handed neutrino 
generates lepton asymmetry at the right level to explain the observed baryon asymmetry of 
the universe. A single phase appears in leptonic asymmetry as well as in neutrinoless double 
beta decay, thus providing some hope for testing high scale phenomena via low energy 
experiments. 

The Model 

We work in the context of low energy supersymmmetry, which is motivated by a solution 
to the gauge hierarchy problem as well as by the observed unification of gauge couplings. 
The gauge group of our model is that of the Standard Model, SU(3)c X SU{2) L x £/(l)y. 
We augment this symmetry with a non-Abelian discrete symmetry A4. This order 12 group 
is the symmetry group of a regular tetrahedron. A<± has a unique feature in describing the 
lepton sector: It has one triplet and three inequivalent singlet representations, thus allowing 
for assigning the left-handed lepton fields to the triplet and the right-handed charged lepton 
fields to the three inequivalent singlets. 

Denoting the three singlets of A4 as (1, 1', 1"), with the 1 being the identity repre- 
sentation, we have V x 1" = 1, V x V = 1", and 1" x 1" = V . Furthermore, 3x3 = 
1 + 1' + 1" + 3 S + 3 a . Specifically, for the product of two triplets we have (ai,a 2 ,a 3 ) x 
(h,b 2 ,h) = (a 1 b 1 + a 2 b 2 + a 3 b 3 ) (1); (ai&i + u 2 a 2 b 2 + u;a 3 b 3 ) (1'); (ai&i + ua 2 b 2 + u 2 a 3 b 3 ) (1"); 
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(a 2 b 3 + a 3 b 2 , a 3 b x + a x b 3 , a x b 2 + a 2 h) (3 S ). Here u = e 2i7r/3 . 

In addition to the A 4 symmetry, we assume a Z 4 x Z 3 discrete symmetry. The Z 4 is 
an i?-symmetry under which the superpotential carries 2 units of charge. The Z 4 and Z 3 
symmetries are broken softly in the superpotential via the lowest dimensional operators. 

The lepton and Higgs fields transform under At x Z 4 x Z 3 as follows. 

L : (3,1,0), e c : (1 + 1' + 1", 3, 0), v c : (3,0,1), £7: (3,1,0), E c : (3, 1, 0), 

# u : (1,1,2), F,: (1,0,0), X : (3,2,0), x':(3,2,l), S 1>2 : (1, 2, 1). (5) 

Here in the fermion sector we have introduced new vector-like iso-singlet fields E and E c 
transforming under the SM gauge group as (1,1,-1) and (1,1,1), respectively, which will 
acquire large masses and decouple. H u and H d are the usual Higgs fields of MSSM, while 
X, x'i &i,2 are all SM singlet fields needed for achieving symmetry breaking. The quark fields 
(Q, u c , d c ) are all singlets of A4 with Z 4 x Z 3 charges of Q(l, 1); w c (0, 0) and d°(l, 2), so that 
the usual quark Yukawa couplings Qd c H d + Qu c H u are allowed in the superpotential. 1 
The superpotential terms relevant for lepton masses consistent with the symmetries is 

W Ynk = M E EiE? + f e UE\H d + ht jk E ie c jXk + \fsrtv c A + \kk^Xk + f v LtfH u .{6) 

Here the flavor structure of the three independent h^ k couplings and the one independent 
fijk coupling can be easily obtained from the A 4 multiplication rules given earlier. 
The Higgs superpotential of the model is given by 

H^Higs = KX1X2X3 + K'six'i + X 2 + X 3 )Si + A^xixaXs 
+ ^s 1± Si + \ S12 S^S 2 + \ S21 S\S 2 + A a22 £f 

+ rfSi + ^ + vAxl + xl + xl)- (7) 

Here x = (XhXz 1X3)1 an d x' = (XiiX^iXs)- The last three terms in the last line of Eq. 
(JJJ) break the Z 4 and Z 3 symmetries softly. The fif 2 terms are the lowest dimensional 
terms that break the Z 3 symmetry softly, while leaving Z 4 unbroken. The (i x term is the 
lowest dimensional term that breaks the Z 4 symmetry softly. Such soft breaking can be 
1 A bare mass term fiH u Hd is not allowed in the superpotential by the symmetries, but the Kahler 
potential, which is assumed to not respect these symmetries, allows a Planck mass suppressed term, C D 
J H u HdZ*d 4 9/Mpi, generating the required n term. 
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understood as spontaneous breaking occurring at a higher scale. We have chosen without 
loss of generality the combination of S± and S 2 that couples to x' as simply Si in Eq. (|7J). 

Minimizing the potential derived from Eq. (JZ|) in the supersymmetric limit, we obtain 
the following vacuum structure: 



(32) 



(3i) = 0; ( X i) = ( X 2) = (X3> = v x , (x' 2 )=V, (xi> = 0; (x' 3 > = 0. (8) 



with v x = —2fi x /\ x , Vg = —fj%/(3\ S22 ), and = (\ S21 lA~ 3A S22 /i^)/(3A S22 A x / s ). Electroweak 
symmetry breaking is achieved in the usual way by (H u ) = v u , (H d ) = f d . We emphasize 
that vanishing of certain VEVs is a stable result, owing to the discrete symmetries present 
in the model. This is important for deriving the MNS matrix of Eq. We observe that 
there are no pseudo-Goldstone modes, as can be seen by directly computing the masses of 
the Higsinos from Eq. (|7j). 

The mass matrices M e E for the charged leptons and M vv c for the neutral leptons resulting 
from Eqs. (JBJ) and (JHJ) are given by (in the notation C = (e, E) M c e (e c , E C ) T ) 







M. 



eE 









h\v x 

h%uiv x 

hoUJ 2 v x 







h%v x 
h\uJ 1 v x 

h\LOV x 



feV d 
f e V d 





M E 










M E 




M,, 









fvVu 

















f»V u 

















fuVu 








fs 2 V s 





fx' v x 


fvV u 








fs 2 V s 








fyV u 


fx' v x' 





fs 2 V s 



feVd 



M E ) 
\ 





h\v x 
h\v x 

\ K^x 




fuV u 



V° 

Since the E and the E c fields acquire large masses, of order the GUT scale, they can 
be readily integrated out. The reduced 3x3 mass matrices for the light charged leptons is 



(9) 
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given by 



M R = Ur 



m P 






m M 



1 

7^ 



111 

1 uj uj 2 
1 w 2 w 



(10) 



where m, = \/3(f e v d v x /M E )hf(l + (hiV x ) 2 ) / M E )~ 1 / 2 . The light neutrino mass matrix is 
found to be 

( 1 
1 - x 2 
x 1 



M light 



m 



V 



11^ 



/ 



(\l + x\ 



where m = f u v u f S2 v s /(f S2 v s - f x >v x ,), and x = -f x >v x >/(f S2 v s ). We define x 
M^ ght can be diagonalized by the transformation M^ s = U*D u lll with 

V2 P*; A 
VI 1 / 
P* is a diagonal phase matrix given by 



\x\e 



iip 



U ^T2 



m 



x 



x 



(12) 



P* 



^{e-* /2 , e -i(^i+^2)/2 ) e -i(fc-Hr)/2^ ^ = arp ^ + ^ ^ = _ x ) _ ^3) 



These Majorana phases will not be relevant for neutrino oscillations, but they will appear in 
neutrinoless double beta deacy and in leptogenesis. The MNS matricx is given by Umns = 
U^U* which has the form given in Eq. (jlj). (In making this identification, we make a field 
redefinition of the left-handed charged lepton fields, = diag.{l, uj, u 2 }e' L .) This is the 
main result of this paper. 



Constraints on neutrino masses 



From Eq. (jTzj) . the expressions for the mass eigenvalues can be inverted to obtain the 
following relations for the parameters |mo|, \x\ and ip = arg(x): 



\m 



\x\ 



cosip 



m 2 



V2m 



im 3 



mlm 2 , + m\m\ — 2m\m\ 



1/2 



— (mt 



m 2 )m 



c 2\p2m\m-i ) [mfm 2 , + m?,m 2 — 2mfm|] 1 ^ 2 
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(14) 



Here m! = |m (l + x)|, m 2 = \m (l — x 2 )\, m 3 = |m (l — x)\. There are restrictions arising 
from the conditions that \x\ be real and | cos^l < 1, which we analyze now. 

Because of the observed hierarchy Am 2 tm ^> Airig, and the requirement of MSW reso- 
nance for solar neutrinos, two possible neutrino mass ordering are allowed, (i) ni\ < m 2 < m 3 
(normal mass ordering) and (ii) m 3 < m x < m 2 (inverted mass ordering). 

If the neutrino masses are strongly hierarchical, m\ <C m 2 <C m 3 , then from Eq. (fT4*|) one 
sees that | cos ^ | < 1 cannot be satisfied, since |cos^| ~ \m 2 j i^2\[2m\)\ ^> 1 in this case. 
Similarly m 3 <C m\ <C m 2 is also not allowed. We find that only two possibilities can arise, 



depending on (a) m\ 



m 2 + m\ (normal ordering), and (b) m 2 — m\ <C m\ + m\ 



(inverted ordering). We consider these cases in turn. 

These conclusions can also be arrived at by analyzing the neutrino mass matrix in the 
flavor basis, i.e., in a basis where the charged lepton mass matrix is diagonal: 



/ 



mo 
3 



3 + 2x - x 2 

2 



V 



-X 

-x 



X 



\ 



x- 



X 



(15) 



O / rp ^ *2 rp 

rp rp^ O rp rp ^ 

When x = — 1 + q with |g| <C 1, we see that entries in the first row and column of Eq. (fTo^ 
become small. This will be the case of normal ordering of masses. 



(a) Normal ordering 

This case is realized when m\ — m\ ~ m\ + m\. The condition | cos^| < 1 can be satisfied 
only if mi/m 2 ~ 1/2. Making expansions in small m\jm 3 and small (2mi/m 2 — 1), we find 



Lr 



1 + 



m\ 
m\ 



cos ip 







2 (2 mi 


- 1 ) 


\ m 2 




m\ I 


2mi 


m 2 \ 


m 2 



- 1 



We see a further restriction that m 1 /m 3 > \2m 1 /m 2 — 1|. 

Neutrinoless double beta decay is sensitive to the effective mass 

\m Q \ 



m, 



x\\l 



x\ 



In the case under study this takes the value 

4 



mp/3 



-mi 



4 / Am| ~ 
6 



1/2 



0.0068 eV 



(16) 



;i7) 



The effective neutrino mass measurable in tritium beta decay m Ue is given by 

1/2 



|2 2 

171, 



(19) 



which in this case takes the value 



m,. 



2 A 2 

3 



1/2 



0.0073 eV. 



(20) 



Here we made use of the fact that ?>m\ ~ AttIq. Finally, the sum of neutrino masses, sensitive 
to cosmological measurements, is given by 



/ Am 2 " 
y^mj ~ 3m! + m 3 ~ 3 I — 



1/2 




+ 


Am »tm 



1/2 



0.064 eV . 



(21) 




Figure 1: Various quantities as functions of \x\ for normal mass hierarchy case, (a) cos^ 
vs. \x\; (b) m\jvn% (dashed) and m 3 /m 2 (solid) vs. (c) m ee (solid (green)), m„ e (dotted 
(red)), m 2 (dashed (black)) and X) 771 ? (dot-dashed (blue)) (in eV unit) vs. (d) RG 

running correction to \U e 3/e\ vs. 
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We have plotted in Fig. la - lc various masses and mass ratios as functions of |x| using 
the exact expressions with central values for Awig and Am 2 atm . These results confirm our 
analytical solutions. 



(b) Inverted mass ordering 

In this case, m\ — m\ <C m\ + m\. The small parameter expansion is different from (a), 
since the denominator of Eq. (16) in cos^ becomes small. Here mi and iri2 are nearly equal, 
and m 3 will turn out to be much smaller than m\. In order to satisfy Am^ <C Am^ m it is 
necessary that cosip ~ \x\/2, which then requires \x\ < 2. Writing 



we find 



cosV = ^(l + g), <2<1, (22) 



A 2 2 2 I 1 2 1 i2/i , nl |2\ 

m = m 2 — m l ~— \mo\ \x\ (1 + 2\x\ )q 



AmL = m l ~ m l - -2|m | 2 |x| 2 • (23) 

It becomes clear that m| < m\ (and thus m| < mf), that is, this case corresponds to an 
inverted hierarchy. With negative q <C 1, the hierarchy in the two oscillation parameters can 
be accommodated. 

In this case the effective mass for double beta decay is given by 

m 3 



[9 + Am r 


1/2 


\ Am* tm ~ 

2 


1/2 


mi 




mi, 





(24) 



with Am^ tm being negative. Here we used the fact that |m | ~ m 3 . The value of m 3 is 
not determined by oscillation data. If m\ ^> |Am^ tm |, we have three-fold degeneracy of 
masses and mpp ~ m 3 . As m\ — > (1/9)| Am^ tm |, the double beta decay amplitude vanishes. 
Furthermore we have for this case 



m„ e ~ mi ~ Y|A7?7| tm | +ml, J2 mi — 2mi + m 3 - 2^/|Am2 tm | + m| + m 3 . (25) 

i 

The three-fold degenerate case is obtained from this case by setting m\ much larger than 
|Am^ tm | (or equivalently, \x\ <C 1). This case also coincides with the leading results of Ref. 

M 

The exact results for various masses and mass ratios are plotted in Fig. 2a - 2c as 
functions of \x\. These results confirm the analytical approximations presented here. 




Figure 2: Various quantities as functions of |x| for inverted mass hierarchy case, (a) cosip vs. 
| x |; (b) mi/m 2 (dashed) and m 3 /m 2 (solid) vs. \x\; (c) m ee (solid (blue)), m Ue ~ m 2 (dashed 
(red)) and X)m» (dot-dashed (black)) (in eV unit) vs. (d) RG running correction to 
|t/ e3 /e| vs. \x\. 

Stability of 

A distinctive feature of the geometric mixing pattern is that = at the scale of A 4 
symmetry breaking, which we have taken to be of order the GUT scale. When running from 
this high scale to low energy scale (Mew), the mixing matrix may change, in particular TJ& 
may not be zero any more. One should ensure that the pattern of Eq. (j3J) is not destabilized, 
which can happen if the induced U e ^ is too large. We demonstrate this stability now. 

The leading flavor-dependent effect of the running from high scale to low scale is given 
by the one-loop RGE (THI 

dM e 1 

This leads to correction, to the leading order, to the entries Mi3 j23 (l — e) and M 33 (l — 2e) 
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with e ~ Y T 2 \n(M G u T /M EW )/32TT 2 . One obtains to order e, 



\ex 




\x 


+ cos ip + % sin ip\ 


3v/2 


cosil>(\x 


+ 2 cos V) | 



(27) 



One obtains the induced \U e z\ for the normal and inverted hierarchies by inserting the cor- 
responding expressions for cosip and \x\ given earlier. 

The results are shown in Figs Id (normal mass ordering) and in 2d (inverted ordering) 
where we plot \U e 3/e\ as a function of \x\. We see that the induced \U e s\ is small, too small 
to be measured by near future experiments for the normal mass hierarchy case in the whole 
allowed \x\ range. For the inverted mass hierarchy case for \x\ larger than about 0.2, [U^l 
remains small. For smaller values of \x\, with e of order one (corresponding to Y T « 1), 
\U e s\ can be as large as 0.1 which may be measured in the future. In this case, all three 
neutrinos are nearly degenerate and the cosmological mass limit on neutrinos will be nearly 
saturated. We conclude that the structure of the mixing matrix derived is not upset by 
radiative corrections. 



Leptogenesis 



Leptogenesis occurs in a simple way in this model via the decay of the right-handed neutrinos 
|14j . The heavy Majorana mass matrix of v c is given in the model as (see Eq. (jOJ) 

" 1 -x- 

M u c = M R 10. (28) 
.-x 1 . 

The Dirac neutrino Yukawa coupling matrix is proportional to an identity matrix at the 
scale of A 4 symmetry breaking which we take to be near the GUT scale. The v c fields 
will remain light below that scale, down to the scale Mr. Renormalization group effects in 
the momentum range Mr < fi < Mqut will induce non-universal corrections to the Dirac 
neutrino Yukawa coupling matrix. Without such non-universality no lepton asymmetry will 
be induced in the decay of right-handed neutrinos. The effective theory in this momentum 
range is the MSSM with the v c fields. 

From the renormalization group equation 



dt 



1 

167P 



(29) 
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where W = eY e LH d + u c Y u LH u + .., we obtain at the scale M R , Y v = x diag(l, 1,1-5) 
with 5 ~ (F r 2 /167r 2 )ln(M GUT /M/j). Y® is the value of the universal Dirac Yukawa coupling 
at the GUT scale. 

We diagonalize M yC by the rotation v c = OQN, where 

/ 1 lx 

, Q = diag{e' l(l)l/2 , 1, e"^ 2/2 } (30) 



y/2 



\/2 
V-l 1/ 

so that the N fields are the mass eigenstates with real and positive mass eigenvalues: M N = 
M R x diag(\l + x\, 1, |1 -z|). 

In the basis where the heavy v c fields have been diagonalized, the Dirac neutrino Yukawa 
coupling matrix takes the form Y v = QO T Y u , so that 

1 + (1 - 5f - e ^-<M/ 2 {(l - 5)2 _ i} 

2 



vyt- l-^°| 2 



L — e 



i(</>i — 02 



)/ 2 {(l-5) 2 - 1} 



1 + (1 - 5) 1 



(31) 



The CP asymmetry arising from the decay of the field N { is given by 

-1 1 



8vr [y,yj 



where 



/(I/) = 



+ ^00 



1 + y 



(32) 



(33) 



,y - 1 y 

As y » 1, f(y) - 3/^y. 

In the normal hierarchy case, mi/m 3 = Mi/M 3 , so the lightest iV field is iVi. In this case 
we have 

-3IK 012 



ei = 



8tt 

±3|Y,°P 2 /mi 



5 2 ( — ) sinU 2 

m 3 , 



8vr 



m 3 



1 - 



(2m 1 /m 2 - 1) 
(m!/m 3 ) 2 



1/2 



(34) 



To see the numerical value of e±, we note that |Y°| can be of order one, 5 ~ (0.1F T 2 ), 
and mi/m 3 ~ [Am 2 olar /3Am 2 tm ] 1 / 2 ~ 0.1. For very large value of tan/3, Y r ~ 1, and 
we find ei ~ 10 -4 . Even for moderate values of tan/3 ~ 20, we find that ei ~ 10~ 6 
is possible. The negative sign will also ensure the correct sign of baryon asymmetry. The 
induced lepton asymmetry is converted to baryon asymmetry through electroweak sphaleron 
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processes. The baryon asymmetry is given by Yg ~ —Y^/2, where Yi = KE\/g*, where 
g* ~ 200 is the effective number of degrees of freedom in equilibrium during leptogenesis, and 
k is the efficiency factor obtained by solving the Boltzman's equations. A simple approximate 
formula for k is jTHj 

(35, 



K ~ 10" 

where 



mi eV 



v 2 



mi = ^Mii (36) 

For 8 ~ 0.1 and Mi ~ 10 14 GeV, we obtain Yg ~ 7 x 10~ n , in good agreement with 
observations. 

For the case of inverted mass hierarchy, A3 is lighter than N±, so we focus on 63. It is 
given by 



4tt Vmi/ ,^l + 2|x| 2 
Again we see that reasonable lepton asymmetry is generated. 

In summary, we have presented a class of renormalizable gauge models based on the 
non-Abelian discrete symmetry A4 which realize the geometric neutrino mixing pattern of 
Eq. (@J) naturally. The resulting constraints on the neutrino masses have been outlined. 
We have also highlighted an intriguing connection between high scale leptogenesis and low 
energy neutrino experiments. 
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